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Abstract 

The relative entropy of the massive free bosonic field theory is studied on 
various compact Riemann surfaces as a universal quantity with physical signifi- 
cance, in particular, for gravitational phenomena. The exact expression for the 
sphere is obtained, as well as its asymptotic series for large mass and its Tay- 
lor series for small mass. One can also derive exact expressions for the torus 
but not for higher genus. However, the asymptotic behaviour for large mass 
can always be established — up to a constant — with heat-kernel methods. It con- 
sists of an asymptotic series determined only by the curvature, hence common 
for homogeneous surfaces of genus higher than one, and exponentially vanishing 
corrections whose form is determined by the concrete topology. The coefficient 
of the logarithmic term in this series gives the conformal anomaly. 
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1 Introduction 



The entropy of a statistical model relative to its critical point has been shown to be an 
interesting quantity in field theory, especially in regard to the renormalization group 
JT|]. On the one hand, it exhibits better behaviour than the free energy when the ultra- 
violet cutoff is sent to infinity and, on the other hand, it is monotonic with the coupling 
constants, unlike the free energy This second property makes it suitable to embody 
the irreversible nature of the renomalization group, which can in particular be sub- 
stantiated in a finite geometry as monotonicity with respect to its characteristic scale 
0. The computation of the relative entropy for various models on a cylinder clearly 
shows its monotonicity M. The cylinder is appropriate to illustrate finite size-effects 
but it may not be the finite geometry of choice in the context of renormalization-group 
irreversibility There is a celebrated result on renormalization-group irreversibility 
in two-dimensional (2d) field theories, Zamolodchikov C theorem. The monotonicity 
theorem for the relative entropy on the cylinder, once it is conveniently formulated, 
resembles Zamolodchikov C theorem . However, this resemblance can hardly lead to 
a direct relationship, since the proof of Zamolodchikov C theorem demands rotation as 
well as translation symmetry. In other words, that proof demands to consider a max- 
imally symmetric space, namely, the sphere, the plane or the hyperbolic plane. Both 
the sphere and the hyperbolic plane possess an infrared scale, the curvature radius, 
but only the sphere is finite and therefore the computation of the relative entropy on 
the sphere is of particular value. 

From a different point of view, the calculation of the spectrum of the Laplacian 
operator on general compact Riemann surfaces has held the interest of first mathe- 
maticians || U and second physicist for some time. The partition function of the 
bosonic massive field theory on a compact Riemann surface is a global object which 
can be constructed from the knowledge of the spectrum of the Laplacian operator. 
Since this spectrum is highly dependent on the topologic and geometric properties of 
the Riemann surface, these properties are reflected by the partition function. However, 
this function, or say the free energy, is ultraviolet (UV) divergent and hence ill defined. 
Fortunately, the relative entropy of the 2d bosonic massive field theory is UV finite 
and so it is likely to have a role in the geometrical characterization of a Riemann 
surface. This characterization consists of local parameters, namely, the curvature, and 
global parameters, specifying the boundary conditions and topologically significant. 
The study of homogeneous surfaces will provide insight into the dependence of the 
relative entropy on both types of parameters. 

Generically, the relative entropy is not related to the quantum field theory entropy, 
but there is a direct relation on the torus (or cylinder) geometry. The relative entropy 
is a geometric entropy, of the sort already considered in connection with the entropy of 
black holes 0. Indeed, the geometry relevant to this case is that of the cone, which is 
non-compact and is actually related to the cylinder geometry, and hence to the usual 
quantum field theory entropy, as were analyzed in Ref. 0. One can expect that the 
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relative entropy for homogeneous spaces will be applicable in a cosmo logical context, 
once suitably generalized. Thus the results of this paper must have some bearing on 
entropic considerations in de Sitter and anti-de Sitter space-time. An attempt at in- 
troducing the maximum entropy principle in quantum cosmology has been made in 
Ref. |§. On the other hand, the application of scaling and renormalization group con- 



cepts in gravitation M and cosmology [10| is gaining momentum. Therefore, it seems 



interesting to study of properties of the entropy relative to the scales defining some 
curved space. Furthermore, the role of the relative entropy as a function monotonic 
with the renormalization group may as well have some relevance in modern theories of 



quantum gravity, as recent work seems to indicate JTT . 

Therefore, our main concern here will be the computation of the universal relative 
entropy of the 2d bosonic massive field theory on homogeneous and compact Riemann 
surfaces, for its own sake and with a view to its application in connection with Zamolod- 
chikov C theorem. For the plane and cylinder, the relative entropy has been computed 
in Ref. 0. Since we are now concerned with compact Riemann surfaces, we shall first 
focus on the simplest case, namely, the sphere. In general, the topological classification 
of compact Riemann surfaces is given by their genus, that is, the number of handles in 
a three-dimensional embedding. The genus of the sphere is zero, of course. Compact 
Riemann surfaces of higher genus, with zero or negative constant curvature, will also 
be considered here, even though they are not globally isotropic, since they are derived 
from the plane or the hyperbolic space by imposing boundary conditions on a finite 
domain which break rotation invariance. The case of zero curvature is the torus and 
is actually related to the cylinder, treated in Ref. @. Compact Riemann surfaces of 
genus g > 1 are always related to the hyperbolic plane. The spectrum of the Laplacian 
on them is too complicated to allow derivation of closed expressions for the relative 
entropy. Thus, the case g > 1 will be discussed summarily and only some general 
properties of the free energy and the relative entropy will be extracted. The method 
used to obtain these properties, namely, the heat-kernel method, is however of general 
interest and we shall dedicate considerable attention to it. 

A large literature has been devoted to the computation of vacuum energy densities 
on various manifolds, mostly in regard to field theory in curved space-time and to the 
Casimir effect. This energy is divergent, of course, and needs regularization. The usual 



technique is the zeta-function regularization p2| , [13| , [14| , as introduced earlier in the 
mathematical literature @]. It is related to the heat-kernel representation that will 
be utilized here by a Mellin transform. We will see how this relationship materializes 
for g > 1 Riemann surfaces in the last section. [] However, since the relative entropy 
is a universal quantity we do not need to bother with prescribing any regularization 
method and we shall only do it to make connections, for instance, with the conformal 
anomaly or with partition functions in string theory. 



1 For a comprehensive review of all these techniques, see Ref. |Tj 
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2 The relative entropy of the sphere 



The free energy W, namely, minus the logarithm of the partition function, for the 
cutoff bosonic massive field theory can be expressed as [|l|, [| 

W[m, A] = — In Z[m, A] = - J> P —^. (1) 

v 

The set of momenta to be summed depends on the type of geometry and is such that 
p 2 < A 2 . If we try to remove the cutoff we see that W is UV divergent, and hence it is 
non universal. The eigenvalues of the Laplacian on the sphere are well known, namely, 
p 2 = I (I + 1)/R 2 , where R is the sphere radius. Thus, the sum over momenta is a sum 
over I and we write 

W[r, A] = 1^(21 + 1) l n '" ( ^' )2 +r . P) 

where we have introduced the dimensionless coupling r = (mR) 2 and the UV cutoff is 
related to the maximum value of I, A w l m &x/R- We have also removed the zero mode 
I = from the sum, which is not allowed for r = 0, and we have taken the degeneracy 
into account by the factor 21 + 1. 

From the previous expression of W one can obtain the relative entropy as a universal 
quantity, that is, as a convergent series in the limit Z max — > oo. However, the sum of 
that series is very hard to carry out, so we choose another way. One can lower the 
degree of divergence of W by taking derivatives with respect to r. In fact, 

dW 1 ^ 21 + 1 



If 21 + 1 

9 i (i i -n 4-7.' 1 } 



dr 2^ l{l + l)+r 

which is still logarithmically divergent. We can remove this divergence just by sub- 
stracting its value at r = and write, in the limit / max — > oo, 



dW dW 



dr dr 



= - T -f +±1 . (4) 



Hence, we define the function U{r) := (dW/dr) — (dW/d^^o, which turns out to be 
computable in terms of the digamma function if>(x). The full expression is rather long; 
it is given in an appendix. 

The function U(r) has interest on its own, since it is the substracted energy, but 
we are more interested in the relative entropy. This can be expressed in terms of U(r) 
as 

S(r) = W(r) - W(0) - r ^YM. = f m s ) ds _ r xj( r \ ( 5 ) 

dr J 
o 
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Unfortunately, this integral cannot be done in closed form. However, it is possible to 
establish the behaviour of S(r) for small or large r. For large r the correlation length is 
much smaller than the radius of the sphere and the result in the plane, S(r) = r/(8ir), 
should be relevant. As R — > oo the sum over I can be substituted by an integral, 



oo 

oo „ 



d 2 p 

W) 2 



where the area of the sphere is A = 4ttR 2 . We see that to compare with the value in the 
plane we must multiply this value by Air. To extract the dominant large r behaviour 
of U(r), we use i/j(x) « lnx. A lengthy but straightforward calculation shows that 
U(r) ~ — lnr/2. Hence, using first an integration by parts, 



r 

S{r) = -J sU\s)ds « T - 



as expected. 

The full asymptotic expansion of S(r) near infinity results from that of U(r), which 
in turn can be worked out with help of the known asymptotic expansion of ip{x). 
However, this procedure yields the asymptotic expansion of S(r) only up to a constant, 
because the condition S(0) = cannot be implemented. It begins as 

S{r) ~ - 73 - tt-—^ + O - . 6 

2 3 lor 105 r z 315 r 6 r 

The presence of a subleading term In r was to be expected, owing to the existence of a 
conformal anomaly for any Riemann surface, coming from the logarithmic divergence 
of the critical free energy. The logarithmic term of the free energy turns out to be 

W = -^HA 2 R 2 ), (7) 

where A is the UV cutoff, R the size, and x is the integral of the curvature divided 
by 47r, equal to the Euler-Poincare number of the Riemann surface, according to the 
Gauss-Bonnet theorem. The calculation of the conformal anomaly from the logarithmic 
divergence of the critical free energy is probably very old but was popularized by the 
development of string theory ||. It was further discussed in Ref. ||13|| , in the context of 



the Casimir energy. In the presence of mass, a term proportional to ln(ra /A ) must 
appear when the correlation length m _1 becomes smaller than R. Hence, it is not 
surprising to have the term lnr in the previous expansion. The concrete way in which 
it appears will be explained in section 4, when we consider the heat-kernel derivation 
of the asymptotic expansion. This is a much more effective method to find the large r 
behaviour, capable of providing the generic form of the coefficients of the asymptotic 
series for arbitrary Riemann surfaces. 
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A plot of S(r) is shown in Fig. 1, in comparison with the asymptotic behaviour 
given by just the two growing terms in the previous formula (||). The agrement is quite 
remarkable, even almost down to r = 1. The numerical value of the missing constant 
is approximately 0.254381. 




2 4 6 8 10 0.20.40.60.8 1 



Fig. 1: S(r) on the sphere compared to the first two terms of 
the asymptotic formula (^) and to the Taylor expansion (|8]) 

The small r behaviour of S(r) is given by the power series expansion near r = 0, which 
is easy to derive since S'(r) = — r U'{r). We obtain that 

S( r ) = \ r 2 + I [2 + ^(1)] r 3 + A [12 + 5^(1) - ^(2)] r 4 + O(r) 5 , (8) 
4 6 lb 

where ip^ n \x) is the polygamma function. The behaviour provided by the series trun- 
cated to this order is compared with the total S(r) in the second plot of Fig. 1. The 
radius of convergence of the Taylor series is determined by the singular points of U(r). 
Since the only singularities of ip(x) are simple poles for non-positive integers, it is easy 
to see that the singularity closest to r = is a simple pole at r = —2, and hence the 
radius of convergence is 2. 

3 The relative entropy of the torus 

We consider the torus as a rectangular box with periodic boundary conditions, which 
is the natural finite geometry in many applications. One can slightly generalize the 
boundary conditions by considering periodicity along two non-orthogonal directions, 
that is, by letting the box be a parallelogram. Although a parallelogram is in principle 
equivalent to a rectangle by an affine transformation, this is only true in real geometry, 
because in complex geometry such transformation is not allowed. Nevertheless, we shall 
consider a rectangle for simplicity and deduce the more general form from holomorphic 
factorization. 
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The partition function on the torus in the critical theory, m = 0, is essentially 
the modulus of Dedekind's function r/(r). The classical proof of this result involves 
the use of the proper time representation and Poisson resummation, after analytical 
continuation in the manner of ^-function regularization, since the partition function 
is UV divergent. This method can be extended to the non-critical theory [16| (also 
fl7fl ). However, we favour a method similar to the one used for the sphere, where we 
calculate the energy U(r). The substraction of the critical value will not be necessary 
since the UV divergence lends itself to straightforward identification. 

Let us L and M denote the periods in the horizontal and vertical directions, re- 
spectively. Then (see appendix) 

oo oo 

W(r) = Y E <0 + E ln t 1 - e " Me( °] + C ' ( 9 ) 

!=— oo £=— oo 



where we have introduced the one-boson energies e(l) = \J (2n l/L) 2 + m 2 and where 
C is a constant irrelevant for the relative entropy. 

Now, we may notice that the previous expression for W can be interpreted as the free 
energy of quantum Id bosons confined in a segment of length L at finite temperature 
T = 1/M and constitutes a slight generalization of the expression for the cylinder 
considered before [0. To prevent a divergence when r = we remove the zero mode 
I = and write 

oo oo 

W(r) = M E e(0 + 2 E ln t 1 - e~ M < l) ] . (10) 
i=i i=i 

The UV divergences concentrate in the Id vacuum energy Eq = X^i e (0- However, 
for r = they can be (^-regularized to give -Eo = (2tt/L) £(— 1) = — 7r/(6 L). The other 
term is the free energy of the bosonic excitations of the vacuum and is finite. For 
r = it combines with E to yield W = — ln Z = ln 77(g) 2 , where rj(q) is the Dedekind 
function, q = exp2iriT and the modular parameter is r = i (M/L). In the case of a 
parallelogram the modular parameter is, of course, complex and Z = 1/ [77(g) 77(g)]. 

Most of the discussion about the relative entropy of the cylinder in Ref. holds 
as well for the torus. Hence, we change the notation for the vertical period, M — > [3, 
in accord with the Id thermodynamic interpretation. The specific relative entropy is 
related to the quantum Id entropy of free bosons, S q , as 

S(r;L, P ) = S { r;L)-^ + ^, (11) 

where 

or r \ / n fn t\ de o( r '> L ) ( j \ E (r;L) 

S{r;L) = e {r;L)-e {0;L)-r — , e (r; L) := . (12) 

The quantity e (0;L) = -tt/(6L 2 ) is dual to -vr/(6/3 2 ). Note that 

S(r;L) = lim S(r;L,0) 

0— >oo 
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is the relative entropy of the cylinder. Of course, lim^oo S(r; L) = Tir/2. While 
S(r;L,P) is modular invariant, the Id quantum entropy S q comes only from the free 
energy of the bosonic excitations, namely, the second term in (|10|) , and is not modular 
invariant. S(r;L,/3) can be computed by the heat-kernel method, studied in next 
section, in terms of a double series of Bessel functions. Its cylinder limit coincides 
with the single series of Bessel functions computed in Ref . [[| . On the other hand, it 
is feasible to obtain the perturbative expansion of W near m = [|l(J], and from it 
the expansion of S(r; L, j3). The asymptotic expansion for large m is just given by the 
value of S for the plane and it comes from the first term, S(r;L), in (|TT|), since the 
second one decays exponentially with L. A more precise description of this asymptotic 
behaviour is provided by the heat-kernel method. 

The quantum Id entropy of free bosons, S q , can be related to the entanglement 
entropy of black holes in 1 + 1 dimensions. The Euclidean version of Rindler space 
is just the punctured plane, which is conformally diffeomorphic to the cylinder. In 
addition, the entanglement entropy of this space can be computed with a path integral 
that reduces it to an ordinary entropy, when considered on the cylinder [[/]]. Therefore, 
one can use the formulas above, taking into account that the conformal transformation 
implies that the modular parameter becomes a function of the quotient of the two 
cutoffs for the radius, namely, i/r = L/f3 = ln(ri/r ). The complete expression of the 
cylinder entropy can be deduced from the results in Ref. . The first term of the low 
mass expansion agrees with the one calculated by other methods in Ref. M . 



4 Heat-kernel techniques 

The computation of thermodynamic quantities for the sphere or the torus, regardless 
of ulterior difficulties, begins with the preliminary step of determining the eigenvalues 
of the Laplacian. Unfortunately, even this preliminary step cannot be taken for an 
arbitrary compact Riemman surface. Therefore, one cannot start from a more or less 
formal but explicit expression of those quantities. Nevertheless, the situation is not 
hopeless for there is a powerful method to extract information on thermodynamic 
quantities, namely, the heat kernel representation of Green functions || £|, |[ ||. We 
now study this method, for the purpose of applying it to g > 1 Riemman surfaces and 
for its own sake: actually, it provides useful insight for the sphere and the torus. 

First, let us recall that 

W(m) = ~Tr 3S lnG(x,y), U(m) = ^Tr x ^G(x,y), (13) 

with G(x, y) the Green function of the Helmholtz equation, namely, (A + m 2 ) G(x, y) = 
8(x, y). If we introduce the Green function of the heat equation 

A,= g, (U) 
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called the heat kernel, K(x,y;t), we can express G as a Laplace transform, 

oo 

G r (x,y) = j dte~ rt K(x,y;t), (15) 
o 

after defining r := m 2 . Furthermore, 

oo 

U(r) = ^J dte- rt Tr^K{x,y-t), (16) 
o 

and integrating over r, 

oo 

W(r) = - l -l je- rt Tr^K(x,y;t). (17) 



o 



In the plane 



K(x, y;t) = — exp 



1 \x — y 



2 



4vrf lJ At ' 

so the specific energy is 

oo 

"M-W*.^., (19) 



A 2 J Aut 
o 

with A = j d 2 x, the total area. To avoid the logarithmic UV divergence at t = we 
may take a further derivative, 



>3C 



™ ±f*e-"= * (20) 



A 87r / 8vrr 



o 



After integrating twice over r, 

W(r) _ 1_ 
A ~ 



(r lnr + C 1 r + C 2 ), (21) 



where Ci and C*2 are two integration constants that are actually infinite due to UV 
divergences. Of course, this expression agrees with the one derived in 0] using UV 
regularization in momentum space. 

In a general curved surface Trg^K(x, y;t) — f d 2 x K(x, x; t) and it is the integrated 
heat kernel for coincident points. Note that K(x, x; t) explicitly depends on the point 



x. However, for homogeneous surfaces it becomes independent of it. According to ([16]) 



^T = \l dte- rt K(t), (22) 



o 
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where K(t) is the heat kernel for coincident points. If the eigenvalues of the Laplacian 
are available, A(fi n = jnipn, the heat kernel is just 



oc 



K(t)=J2 e7nt - ( 23 ) 



n=0 

The first eigenvalue is 70 = 0, corresponding to the constant solution; the others are 
negative and ordered by their absolute value — in fact, 7! < —1/4 M. Substituting for 
K(t) in the integral representation fl22|) we recover the expression for the energy used 
before, 

^ = -E— ■ (24) 

n=0 

For example, if we consider the torus, 7„ = — (27r) 2 |/c| 2 , where k G Q*, the lattice 
dual to Q, the one definig periodicity on the torus. Then we can use the Jacobi identity 
to write 

K(t)= ^exp[-( 27 r) 2 |fc| 2 t] = A^ exp[ _tf]. (25) 

ken* wen 

The latter sum is 1 + 0[exp(— l/t)], reproducing the result for U in the plane plus 
corrections vanishing exponentially when r —>■ 00. 

Generically, the heat kernel K(t) admits an asymptotic power series expansion at 
i = 0@, 

+ (26) 

in terms of some numerical coefficients determined by the Riemann curvature. This 
expansion, upon integration over t, leads to asymptotic expansions of thermodynamic 
quantities for large r. For the hyperbolic plane, it can be calculated from the integral 



of exp — v t times the density of states of frequency v |T9], |20 



2tt 

It yields 



00 

K{t) = / dvv tanh(7rz/)e^ *. (27) 



t/4 [ °° 1 _ 2 l-2n ] 

m = ^-^~^-B^, (28) 



where Bm are the Bernoulli numbers. The expansion for the sphere is related to it by 
analytic continuation from the hyperbolic to the elliptic geometry and has an additional 
(— l) n in the sum. 
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From the expression (|22|) for the energy we can derive the heat kernel integral for 
the specific relative entropy, 



S(r) 
A 



oo 

~ / j[l-(l + rt)e- rt ]K(t), 



(29) 



which can be easily shown to be convergent at t = 0, taking into account the behaviour 
of the heat kernel for small t (0). For a compact surface, it also converges as t — > oo if 
we remove the zero mode, since then K{t) decays exponentially, according to Eq. (^3|). 
To obtain the asymptotic r —>■ oo expansion we could be tempted to just substitute 
the aymptotic expansion ( p6f ) into this integral. However, we would encounter that the 
coefficients are given by integrals divergent as t — > oo. The problem is that we cannot 
prolong the aymptotic expansion (|6|) to t = oo. However, the aymptotic expansion 
(P^) can be substituted into the integral for S"(r), 

oo 



A A 2 

b 

to provide its asymptotic expansion r — > oo, owing to Watson's lemma [^T|. This 
expansion can in turn be subjected to indefinite integration to yield S(r), of course, 
up to a constant: 

S(r) „ 1 ( , ^ (n + 1)! A „ n 

We can obtain the whole asymptotic series of S(r) for the sphere by just taking the 
coefficients from the corresponding asymptotic series of K(t), given by Eq. fl28|) with 
(— l) n inserted. The coefficients of its succesive terms coincide with the ones in the 
previous expression (|6|) except for the one of the logarithmic term. This term, which 
was already remarked, is particularly interesting and we now discuss it in more detail. 



4.1 The logarithmic term in the asymptotic expansion of S(r) 
and the conformal anomaly 

The logarithmic term in the asymptotic expansion of S(r) comes from the constant 
term in the asymptotic power series of K(t) (f26[), namely, ai/(4ir). In particular, 
formula fl28|) yields a\ = ±1/3, for the sphere and hyperbolic plane, respectively. If we 
substitute this constant for K{t) in fl30D we obtain 

oo 

A ~ A ~ 2 J Un~ Svrr' W 

o 
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(Notice the slight abuse of notation, for we deal here with truncated quantities.) For 
a curved compact homogeneous surface the Gauss-Bonnet theorem implies that A = 
2vr |x| = 47r |<7 — 1|, where \ is the Euler-Poincare number. Then the coefficient of the 
logarithmic term of the total value of S(r) is 

gi A = X_ = 1 ~ 9 
8vr ~ 12 ~ 6 ' 

only related to the topology of the Riemann surface. This result also holds for compact 
surfaces of variable curvature. 

We can obtain as well the coefficient of the logarithmic term of W. Upon integrating 
U'(r) twice over r, we conclude that 

EM^^ + Cr + ft. (33) 

A 07T 

Hence we can try to connect with the critical value stated in (|7|) as follows. Let us 
return to the heat-kernel representation and perform the divergent integral for W(r), 
Eq. (|lTD, with a lower cutoff: 



W{r) ~ ^ j d {e-^ = -^m m ^) = -^[-H^e)- 7 + 0( m \)}, (34) 



A 8tt 



where T(0,x) is the incomplete gamma function. This agrees with the previous result. 
We may interpret the cutoff in t as an UV cutoff, e ~ A" 2 . Reinstating the radius of 
curvature we can then split m 2 e so that 



ln(m^e) = lnr - \n(A z IT). 

The second term agrees with the critical value, Eq. (|7|), but the first term diverges as 
r — ► 0. In fact, this limit is not meaningful since the logarithmic term arises in the 
asymptotic expansion for large r. Neither is it meaningful to utilize the asymptotic 
expansion of K(t), Eq. (p6|), to obtain the small-r behaviour: for r = the integrals 
for thermodynamic quantities are determined by the entire range of t and not just by 
its asymptotic behaviour, whether it be for small or large t. Nevertheless, Eq. ( p6|) can 
still be utilized to determine the form of the divergence of W(0) at t — 0. Indeed, 

oo 

A f Hi 

W[0,A] = -- J jK(t) (35) 

A- 2 

implies that 

A 9W t A] = -AK{A- 2 ) = - A [A 2 + ox + 0(A- 2 )] . (36) 

OA 47T L J 



12 



In the infinite-cutoff limit, AR 3> 1, we can discard the negative powers of A and this 
equation is none but the conformal anomaly equation for a rigid scale transformation 
(recall that a\ is just the Riemann curvature divided by 3). Integration over A yields 
the desired term, 

4 " 1 \n(AR) = -§ \n(AR), 



4vr 6 
on a compact surface. 

Notice however that W(0) can also be calculated with zeta-function regularization 
and that it yields a finite result || . This may seem puzzling, for we have then lost track 
of the conformal anomaly. In any event, W(0) is a pure number and its value is to some 
extent irrelevant. Nevertheless, in the massive case, zeta-function regularization also 
yields a finite W(r), as we shall remark in the next section. But then the conformal 
anomaly can be extracted from the asymptotics of W(r) for large r. Moreover, with 
any type of regularization, the asymptotics of S(r) for large r provides the conformal 
anomaly. This conclusion is not particularly interesting for the massive bosonic field 
theory, the only one considered here, but it may be very interesting for interacting 
theories. 

We have concluded that the coefficient of the logarithmic term in the entropy must 
be the opposite of the conformal anomaly, that is, x/12. While we had obtained —1/3 
in Eq. (|6]) now we get xjVl = 1/6. This discrepancy stems from having suppressed 
the zero mode in the calculation of U(r) for the sphere. This contributes 1/2 in any 
compact surface and 1/2 — 1/3 = 1/6. We may make a little disgression here and 
recall that the calculation of the small-t behaviour of tensor Laplacians on compact 
Riemann surfaces provides a proof of the Riemann- Roch theorem || . In the case 
of the scalar Laplacian, A = <r 8, this theorem states that 

1(d) = dimKer<9 — dimKer^ = 1 — g, (37) 

where the adjoint tr = V* is the covariant derivative on forms. Then KerV* is the 
space of Abelian differentials, with dimension g and Ker B is the space of holomorphic 
functions, that is, constants, with dimension 1 (the zero mode). Their difference is 
1 — g — x/2. According to our discussion on the presence of logarithmic terms in 
r, we can interpret the zero mode as the logarithmic term for r — > whereas the 
logarithmic terms found in the asymptotics r — > oo are related to the existence of 
non-trivial boundary conditions and hence with the existence of Abelian differentials. 
They substract from 1 the right number g such that the difference is proportional to 
the Euler-Poincare number. Indeed, we observe that, e.g., for the torus the zero-mode 
term In [l — e~ Lm ~\ ~ ln(m L) as m — * 0, while it vanishes exponentially in the opposite 
limit, m — > oo. This is in accord with the torus being flat, so that \ = 0. 
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5 Thermodynamic quantities for compact Riemann 
surfaces of higher genus 



A compact Riemann surface of g > 1 can be characterized by its fundamental group. 
When this surface is represented in its covering space, the hyperbolic plane, it gives rise 
to a tesselation of it, whose tiles are identified by a discrete group of motions isomorphic 
to the fundamental group. Since the total group of motions of the hyperbolic plane 
is SX(2, H), that group is one of its discrete non- commutative subgroups, which are 
called Fuchsian groups. This construction is analogous to the construction of the torus 
by factoring the plane by a lattice Q, where the fundamental group is Z x H. We have 
seen that there is an expression for the heat kernel of the torus as the kernel in the 
plane plus a series of powers of exp(— 1/t) (|25f) . There is a non-commutative analogue 
for g > 1, namely, the Selberg trace formula 0J], 

oo 

m = -J ^^.-^SS exp h ^f ]. 

o n=1 {7} 11 

(38) 

In this formula {7} are the primitive conjugacy classes of the Fuchsian group and Z 7 is 
the length of the shortest geodesic along the corresponding non-contractible loop. The 
first term is just the integral representation of the heat kernel in the hyperbolic plane 
H fl2"7p. It admits a different representation, more useful for computations, derived 
from the form of the corresponding Green function Jl], ^j, 

00 

K H (t) = / db———. 

1 sinh (6/2) 







Although the Selberg trace formula is a much more complicated formula than the 
Abelian one Q2"5|), it has the same structure; namely, it is a sum of the part correspond- 
ing to the infinite surface, now the hyperbolic plane, plus 0[exp(— 1/t)] corrections 
due to the boundary conditions. As well as for the torus, these corrections give rise to 
exponentially vanishing terms which do not appear in the asymptotic expansion of the 
entropy around r = 00. Therefore, all the higher genus Riemann surfaces share the 
same asymptotic expansion in r. 



Formula (38) has been widely used to establish the modular dependence of partition 
functions in string theory || |2^]. The partition functions in string theory are those 
of critical theories and Selberg's formula results in generalizations of the Dedekind 
function of the torus to g > 1; they are called Selberg's zeta functions [|J. Although we 
are interested here in the massive case, the procedure to calculate the critical W @, U 
actually applies to the non-critical one as well. It consists of a part corresponding to 
the hyperbolic plane, hence common for any homogeneous surface with g > 1, and a 
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part in terms of Selberg's zeta function, which accounts for the topology of the surface: 

1 



1 H r) = -- C(0) - (4 r + In Z s ( Jr + \ + \ 



(39) 



The zeta-function of the hyperbolic plane is 



Cr(ac) 



p -(r+l/4)i /■ /, p -6 2 /(4t) 







sinh (6/2) 



iJm + 4 r ) s/ " I/2 / * dSJj 6 f+ i 

v 



(40) 



(41) 



Notice that it is a meromorphic function of x with one single pole at x = 1, as long as 
|r| < 1/4. In the critical case the common part can be calculated exactly to yield 



Co(0) = 2G/-1) -ln(2vr) + --4C / (-l) • (42) 
Selberg's zeta functions is defined as 

oo 

Zs(x) = Hl[{l-e-^ 1 ^. (43) 

7 p=0 

The function ( r {x), as a regularization of W on the hyperbolic plane, must lead to 
the asymptotic expansion of the relative entropy provided by Eqs. (ESp and Q2Tf). The 
Selberg zeta function leads to the asymptotically vanishing corrections. 



6 Conclusions 

We have seen that the general structure of the entropy of the free massive bosonic field 
theory on compact homogeneous Riemann surfaces consists of a part corresponding to 
the maximally-symmetric surface, namely, to the sphere, the plane or the hyperbolic 
plane, and a part due to the boundary conditions. The first part can be expressed 
as a complicated function, analytic in r e [0,oo). Furthermore, it has an asymptotic 
expansion around r = oo, which is fully computable. The second part embodies the 
topology and is more delicate to treat, but it vanishes exponentially as r grows and 
therefore does not appear in the asymptotic expansion. The behaviour of the entropy 
for small r — the critical limit — is also calculable, in terms of a convergent series. It 
is strongly dependent on the global parameters defining the topological nature of the 
surface. Indeed, one can observe, for example, that the series for the sphere (§D has 
nothing to do with the one for the cylinder, obtained in 0. The kth term of the small- 
r expansion is easily seen to be proportional to J2 n (—•y n )~ k , beginning with k = 2. 
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But for g > 1 these sums cannot be computed analytically, because the Laplacian 
eigenvalues are not available. 

The monotonic character of the entropy can be checked in our calculations. For large 
m the entropy tends to m 2 /(87r), independently of the topology of the surface, which 
is the value for the plane. However, the topology lets itself be felt in the subleading 
term, proportional to In m, which is actually related to the conformal anomaly. These 
two terms are the only ones divergent as m — > oo. There is also a constant term in the 
asymptotic expansion, which cannot be determined exactly, however. The three terms 
together already provide a good approximation down tomL ~ 1. As we decrease m we 
approach criticality and the entropy becomes very sensitive to large scale peculiarities of 
the surface, as already remarked, but one can check that it always remains monotonic. 

We may wonder how much of the above can be generalized to higher dimensions 
or to interacting field theories. The spectrum of the Laplacian and the heat kernel 



for homogeneous spaces are well known | 20fl . In fact, the heat kernel is simpler in 



odd-dimensional homogeneous spaces than in even-dimensional ones, so the formulas 
for the three-dimensional sphere or hyperbolic space turn out to be simpler as well. As 
regards massive interacting theories, an asymptotic expansion for large m must exist 
in general and, moreover, the leading and subleading terms can also be studied. The 
leading term is always proportional to m 2 for dimensional reasons, and its coefficient 
positive. To determine this coefficient one can use the thermodynamic-Bethe-ansatz 
computation of the universal bulk term of the free energy, as discussed before ||. The 
subleading term is likely to be related to the conformal anomaly, c(x/6) lnm, where 
c is the central charge of the conformal field theory for m = 0. This may provide a 
new way to find the central charge of a model. The small-m behaviour can be studied 
with conformal perturbation theory, see Ref. [0. However, it shall crucially depend on 
the nature of the surface, as well as on the particular field theory. This perturbation 
theory is presumably convergent for strongly relevant perturbations, like that of the 
free bosonic field theory. 
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Appendix 



Substracted energy for the sphere 

We here perform the computation of the substracted energy for the sphere. 



U(r) 



taking into account that 



^[/(/ + l)+r](/ + l) ^[l(l + l)+ r }l 



E 



(44) 



1 

2r~ 



27+ 1 



VI -4r 



. ,3- VT^F, . , 
V>( ^ )+ 1 



E 



vr^ 



3 + y/T=T: 



,(45) 



E 



((l- 7 ) (1 -4r- VI -4r)) 



^ [J (J + 1) + r] (J + l) 
(VT^47 + 2r- 1)^(^4^) ~2r^( 1+v ^ ) 



(-1 + VI -4r) VI -4rr 



we have that 



VT^. 



(46) 



3 + VT^47, 

n o ) 



2 ((1 - 7 ) (l - 4r - VT^47) + (VT^47 + 2r - l) ^( ^v^ ) + 2r ^(±±4^)) 
+ (-1 + VT^47) VI -4r 

The apparent pole at r = must cancel and a careful analysis shows that it does; the 
pole at r = 1/4 cancels as well. Notice that when r > 1/4 the argument of the digamma 
functions becomes complex. Nevertheless, U (r) remains real and is an analytic function 
of r at r = 1/4. 



.(47) 



W(r) for the torus 

Let us L and M denote the periods in the horizontal and vertical directions, respec- 
tively. Then 



W[m] = -\nZ[m] = X - In +( 

l,n=—oo L 



'27m 



+ m 



(48) 
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and 

U(r)-=^= 1 - V I (49) 

lJ ' 2 ^ (l/Lf + (n/M) 2 + r' 1 J 

l,n=— oo 

where r := (m/2ir) 2 . To work out this sum we can use the known expansion of the 
hyperbolic cotangent in simple fractions [23| , 

COth(TTx) = ~ ^—2- ( 50 ) 



n=— oo 



Hence, 

C/(r) = - Y, 7ffifV 2 coth(nMy/(l/L)* + r). (51) 

Now we can obtain W(r) by integration. Notice, however, that the series (|5l|) 
is divergent, so term by term integration is not warranted. However, the series for 
U(r) — U(0) is convergent and one can apply it to it. The ensuing series represents 
W(r) — £7(0) r minus a UV quadratically divergent constant. Of course, U(0)r is the 
UV logarithmic divergence of W(r). Since the UV divergences of W(r) can be easily 
segregated, we can proceed with the integration term by term without further concern: 

Y J /(//L)2 + r coth(nM y/(l/Ly + r) = — J de coth( y e) = lnsmh(y e), 

where we have introduced the one-boson energies e(l) = v/(27r//L) 2 + m 2 . Finally, 

°° M M °° 00 

W(r) = £ lnsinh( y e(l)) = y ^ + £ In [l - e -"*>] + C, (52) 

!=— oo £=— oo £=— oo 

where C = — ln2^]^_ 00 1 is an irrelevant divergent constant. The other divergences 
appear in the first term and are of the form C\ + C2 r, as already remarked. 
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